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LIST OF SYMBOLS 



Im r,) normalized attenuation constant 

C3 ^ 

Re r ) normalized phase constant 

^a ^ 

speed of sound in the air of the pores 

effective speed of sound in the air of the pores (see Appendix) 
speed of compressional waves in the frame structure 
speed of sound in air in free space 
mean diameter of a filament in fibrous structure 

k"/K P-Y, generalized structure factor 

5 d O 

structure factor for air in the pores 
depth of sample of finite size 

instanteous acoustic pressure in the air of the pores 
static acoustic pressure in the air of the pores 
acoustic pressure amplitude in the air of the pores 

acoustic pressure amplitude in the porous frame structure 

to p/$ , generalized acoustic quality factor 
(op+ i $ , generalized particle velocity coupling coefficient 
acoustic particle velocity amplitude in the air of the pores 
acoustic particle velocity amplitude in the porous frame structure 
complex acoustic admittance of air in the pores 

complex acoustic admittance of the porous frame structure 

complex acoustic impedance of the air in the pores 
complex acoustic impedance of the porous frame structure 





LIST OF SYMBOLS 


attenuation constant 

normal incidence intensity absorption coefficient 
phase constant 

instantaneous acoustic density ratio 
1.4, ratio of the specific heats for air 
complex propagation constant 

complex propagation constant for sound waves in the air of the pores 
complex propagation constant for sound waves in the elastic frame 
structure 

bulk modulus for air in the pores 
bulk modulus for elastic frame material 
generalized dilatational coupling coefficient 
apparent density of air in the pores at high frequencies 
effective density of air in the pores at all frequencies (see 
Appendix) 

Po(k$ - l)/fi, effective density in the inertial coupling coefficient 

effective density of the frame structure 

mean density of air in free space 

mean density of the material of the frame structure 

coefficient of kinematic viscosity for air 

specific flow resistance of air in the pores 

dynamic resistance of the air in the pores (see Appendix) 

generalized elastic coupling coefficient 

porosity of the frame structure 

angular frequency of excitation 



I. Introduction 

Highly porous, fibrous materials play an important role in the 
technology of sound control. Such materials are used in a wide variety 
of forms and configurations for special purposes. Some fibrous materials 
such as glass fibers are synthetic and are spun together to form a flexible 
blanket while others appear as open-cell elastic foams. Others consist 
of natural fibers such as wool or hair felt. Most are available in a 
range of porosities and densities. Natural surfaces such as ground cover 
sometimes are modelled analytically, somewhat similar to fibrous bulk 
materials. Frequently, the bulk material is covered by a facing to 
prevent damage, provide mechanical support, and even change appearance. 

Interest in this study is confined to acoustical properties of bulk 
materials and how their physical parameters influence acoustical character! sti 
such as normal incidence impedance and intensity absorption coefficient 
at its surface. No attempt is made to optimize performance in a particular 
application. In spite of the widespread application of both commerical and 
natural fibrous materials and their importance in noise control very little 
basic research has been conducted since the pioneering work of Scott^^^, 
Zwikker and Kosten^^^’ Morse and Bolt^^\ Beranek^^^ and R. W. Morse^^^ 
to mention a few early workers. Their primary motivation was better under- 
standing in applications to architectural acoustics. 

In this report the properties of fibrous materials (blanket, foam, and 

stiff tile) are formulated analytically on the basis of the theory of Zwikker 
( 2 ) 

and Kosten' ' as applied to elastic frames. So, fundamentally, the equations 
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of motion are not new. However, they are recast into a form that identifies 
and hopefully clarifies certain aspects of coupling between sound waves in 
the air of the pore space and compressional elastic waves in the fibrous 
frame structure. Analytical solutions are then obtained for weak coupling, 
which appears to be most valid for highly porous materials, and an attempt 
is made to specify the minimum number of acoustical parameters required to 
characterize the behavior, and to establish any apparent frequency limitations 
in the characterization. The ultimate objective is to develop design procedures 
required in critical applications such as flow duct acoustics. 

The coupled mode theory of Zwikker and Kosten sometimes has been criticized 
as being somewhat ad-hoc and too complicated. While this may be true for 
certain simple applications it is not valid for critical design where specific 
impedance values are needed. One important reason for studying coupled wave 
motion in bulk elastic materials is a better understanding of the behavior 
of duct liner surfaces that may exhibit extended reaction. Flexible fibrous 
materials that exhibit weak coupling are not likely to be locally reacting 
especially if the sound speed in the air of the pore space is much different 
from that of the sound speed of the incident wave. 

Several experimental studies were carried out in an impedance tube to 
establish properties of small samples of a highly porous, fibrous material 
called Kevlar 29. One type had been studied earlier and some of it bulk 
properties were known^^^. Kevlar 29 consist of long, very fine fibers, and 
has excellent sound absorptive properties in bulk form^^\ There are two 
types, one in which the fibers are layered and the other in which the fibers are 
woven to form a more dense and stiffer blanket. These absorptive properties have 
been ascribed^^^ to a very slow phase speed at low frequencies which leads to 
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short wavelengths and more wavelengths per unit depth of sample. The 
attenuation at very low frequencies is about 2Tr nepers per wavelength. In this 
report experimental values of normal incidence impedance and absorption 
coefficient are compared with predictions based on decoupled wave motion. 

The agreement with theory is satisfactory if the so-called structure factor 
(ratio of apparent density of air in the pores to the mean density of air 
in free space) is adjusted to provide a good fit to impedance data. 

A parameteric procedure was established to provide that fit. Future work 
could concentrate on a characterization based on coupled wave motion if low 
frequency impedance data (below 300 Hz) becomes available. The weak coupling 

formulas are derived in this report. 

A comprehensive experimental study had been conducted earlier on the 
acoustical properties of a highly porous, open-cell, flexible, fibrous 
foam called Scottfoam^®^ This material also has excellent absorptive pro- 
perties and a comparable parametric study yielded satisfactory results based 
on decoupled motion. It turned out that an acoustical quality factor index 
for Scottfoam was roughly twice that for Kevlar which makes it a very 
efficient bulk absorber. Microscopic study revealed that the filaments in the cellu 
lar structure have a mean diameter comparable with that of Kevlar, see Table I. 

In the section on parametric studies an initial attempt is made to relate 
some of the material parameters in the theory of Zwikker and Kosten to 
parameters developed by others notably Morse^^^ and Hersh and Walker . Also, 
in the appendix a new derivation of the equivalent density and dynamic 
resistance for porous materials as first proposed by Beranek is presented . 

While not directly relevant to this paper it could have considerable relevance 
if future work on the extended reaction of bulk materialsin duct liners is 
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undertaken. It turns out that these more specialized theories can be 
brought into agreement with the theory of Zwikker and Kosten by suitable 
definitions. 

In the following section a reformulation of the equations of motion 

and resulting characteristic equation based on the thoery of Zwikker and 

Kosten is presented and discussed. The notation is consistent with that used 

(9l 

in a modern textbook such as Morse and Ingard' ' . 

Use of commercial products or names of manufacturers in this report 
does not constitute official endorsement of such products or manufacturers, 
either expressed or implied, by the National Aeronautics and Space 
Administration. 
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II. Analysis 

A. Characteristic Equation 

One of the important questions raised in the introduction is deciding 
whether or not sound waves in the pore space and elastic waves in 
the fibrous frame are coupled in any significant way. Roughly 
speaking one expects maximum coupling to occur when the phase velocity of 
sound waves in the air is equal to that of the longitudinal elastic waves 
in the frame. This coupling is both inertial and dissipational and is caused by 
differences in the wave particle velocities which must be accounted for in 
the balance of forces. Also present is coupling due to different pressure 
forces which act on the volume under consideration and produce changes 
in volume which are expressed in terms of "elastic constants" in the gas 
and in the frame of fibers. This coupling must be accounted for 
in the equations of continuity which conserves the total mass in a given 
volume of material. The elastic constant of the gas is fairly well defined 
but the elastic constant of the frame is not so well defined and hence 

• i • At 

not easily measureable either. Usually a static value is used 
present this matter can only be resolved by doing careful measurements v/ith 
the objective of resolving any discrepancies between measured and predicted 
values of the propagation constants, surface impedance, and acoustical 
absorption properties based on decoupled motion. Therefore, it is worthwhile 
studying wave coupling for the case of flexible, highly porous materials 
consisting of fine fibers and determining what physical properties are 
required to effectively decouple the motion especially in the lower ranges of 
frequencies say below 500 Hz. Since, the air motion is viscous controlled 
over a fairly wide range of frequencies for such materials it is not immediately 
obvious as to exactly what requirements must be satisfied. 
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Thick samples of porous materials are homogeneous for all practical 
purposes and by considering plane wave propagation in one given direction 
the material may be considered isotropic. However, some fibrous materials 
are layered causing the material to be anisotropic. Thus any angular 

dependence of either impedance or absorption coefficient could 
not easily be inferred from either normal incidence measurements or a 
one-dimensional theory. Foam materials on the other hand, appear isotropic 
with no preferred orientation of the intercommunicating filaments and hence are 
probably easier to analyze and project results. 

It was decided to use the basic formulation due to Zwikker and Kosten 
asa starting point in the analysis but to recast the notation in a more 
familiar form and use definitions as presented in a modern text like 
Morse and Ingard^^^. Following linear theory the relationship between the 
changes in air density 6 caused by acoustic pressure disturbance p in the 
pore space is assumed to be of the form 

6 = po<aP 

where po is the mean density of the air and is the bulk modulus of 
the air which, generally speaking, is frequency dependent. For very low 
frequencies and very fine pores it is assumed to take on its isothermal value 
1 /Pq where Pq is the static pressure while at sufficiently high 

frequencies it is assumed adiabatic, i.e. Ka “ |/PoY where y is 
the ratio of specific heats for air. Initially no assumption about the value 
of Ka will be made except that it is treated as a real number. That matter can be 
deferred to a parametric analysis. Section III. 
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The solid frame structure of the porous material is characterized in vacuum 


by a bulk modulus ks and a mean mass density ps in a manner similar 
to that of air. This type of elastic medium supports shear waves but 
only internal coupling effects of compress ional waves is "uiiportant 
for normal incidence predictions. The dependence of Kg on frequency is not 
known, generally speaking, for most fibrous porous materials including foams. 
Initially no assumption about k^ will be made either. 

Another parameter of great importance is the "apparent density" of the 
air in the pores, Pg. Some fibers will move with the air as sound propagates 
thereby increasing its apparent density. Also, the air particle motion is 
constrained by the presence of constrictions and the apparent density is 
increased when air is forced to move through such constrictions. At sufficiently 
high frequencies the density of the air in the pores is assumed to be given 




( 2 ) 


where k$ is called a "structure factor" which supposedly accounts for 
the increase in air density caused by constrictions and blind passages, and 
is the porosity of the bulk material. Note, an increase in 
density caused, by motion of the fibers is not included in Pg . 

That matter is accounted for in a different way, see appendix A. For most fibrous 
materials 1 < p /p < 4.5 . Usually ks is close to unity for fine 
fibered, highly porous materials, the finer the fiber the lower the value 
for a given porosity. Actually kg and must be interrelated but the 
exact functional form is not important here. 

Both internal inertial and frictional coupling can take place between the air 
motion and frame motion as mentioned above. Only internal coupling will be 
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considered even though some wave coupling can take place at both surfaces 

(o) 

if the material is of fipite depth . Boundary value problems will not 

be formulated until specific applications are considered. 

For harmonic compressional waves of the form exp [-i(rx - tot)] the 

important inertial frictional coupling parameter in the theory of Zwikker 
( 2 ) ’ 

and Kosten^ ' is here designated as 



where Pei is a density and $ is a specific flow resistance (rayls per unit 
length). For most fibrous blanket materials ks £ 4. At low frequencies 
where ‘^^01 ^ the internal coupling is primarily frictional. 

Also define an effective density of the structural frame Pg^ as 


es 


= Pc 


n - n) 




( ks - 1) 

Q 


(4) 


where is the density of the bulk fibers. This definition incorporates 

Pg^- where p^^. < <Pg^^ ~ even for porous materials for which 

> 0.95 and k^ ^ 4 . 

The ratio Pq^/Pq is an important parameter in the theory of internal 
wave coupling. It is an index of the potential for appreciable acoustical 
energy storage in the frame as compared to the air in the pores of the frame. 

The flow resistance ^ , Eq. 3, is not independent of either kg or 
If the fibers are layered, a one-dimensional wave model with propagation normal 
to the direction of the fibers suffices, provided $ is specified in the 
same direction. Likewise for propagation along the fibers. No attempt will 
be made here to introduce anisotropy and derive a three-dimensional wave 
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equation. Propagation normal to the fiber corresponds closely to some 
applications. 

The propagation constant r corresponding to one-dimensional wave 
motion.will be complex number 6 + ia whose individual properties are 
important. 

The dynamical equations of motion including the force equation and 
equations of continuity have been reformulated to reflect these definitions 
and also to identify important coupling parameters. (See Appendix B) 

For the elastic frame the force equation can be expressed as 


r Ps = SI [copgg + i$] Vs - [wPg^. + i$] Va . 


(5a) 


By using definitions put forth by Kosten and Janssen^^*^^ the equation 
of continuity becomes 


wKs (1 " Si) Ps = Kii r Vg + <12 r Va 


(5b) 


where 


^11 


s 1 +h 








(1 - Po^a) 


and 

■=12- Si 2 • 

The variables Pg and Pa are the complex acoustic pressure amplitudes 
in the frame structure and in the air, respectively, while and Vg are 
complex particle velocity amplitudes. These variables are taken as averaged 
values over areas large compared to a "mean pore size". The precise way that 
average is computed is not clear. 
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For the air in the pores the acoustic force equation can be expressed 
as 

r Pa = (ojPg + i$) Va - « (wPg^. + i<I>) (6a) 

while the equation of continuity becomes 


Pa = <22 ^ '^a '<^21 ^ 


(6b) 


where <22 = 1 and 

*^21 ■ (^ ‘ ^) " ^o'^a^ * 

Note, that the dilatational coupling parameters tc -|2 and < 2 q simplify 

somewhat if <=. = -^ but 
a ro 

<22 ■ '^12 *^21 "" ^ 

regardless of the thermodynamic process involved. Not further that 


S-ji $22 - S^2 ^21 ^ 


Equations 5 and 6 are expressed more conveniently in matrix form below: 


pressure 



velocity 



0 

■^11 

$12 


p -1 






Ps 

0 

r 

'21 

-S 22 

i nertial _&_f ri ctlonal 


Pa 



^11 

r <12 r 


Vs 

0 

-022 

<21 

r <22 ^ 


Va 

elastic 



dilatational 



= 0 


( 7 ) 
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where in addition the parameters 

(wpgg + i$) 
$22 = wpg + i$ 

a,, = (1 - n) 

11 s 

Ooo = 

cd a 


^*** c / • \ r- 

S21 = Q 

First note that both Eqs. 5 and 6 reduce to their proper forms 
as either ft ^ 0 or ft 1 . Second, the matrix of wave coefficient 

in Eq. 7 has a certain symmetry with the inertial -frictional terms in 
the upper right hand corner, the elastic terms in the lower left hand corner, 
and the dilatational terms in the lower right hand corner. Notably, absent 
in this formulation are any terms corresponding to elastic cross coupling, 

i.e.t ^12 ~ *^21 ~ ^ * 

Observe, however, that the k's do contain elastic constants and the 
relative importance of say <12 depends on both ft and the elastic ratio 
Kg/Kg . It also turns out that = 1 for highly porous, bulk materials 
and both k -|2 and <2] can be quite small compared to unity* in fact 
<12= 0 ii" <a = 1 /Pq ♦ So all in all no significant elastic wave 
coupling is present in the theory of Zwikker and Kosten. Note further that 
since iCn <22 “ ><12 '^21 ~ ^ dilatational coupling that exists is 
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reciprocal. 

There will exist four values for the propagations constant r corresponding 
to the coupled motion and waves travelling in both positive and negative x directions 
Many aspects of the solution for r have been discussed by Zwikker and 
Kosten^^^ and somewhat later by Beranek^^^ and McGrath^ and their 
interpretations will not be repeated here. By setting the determinant of 
the matrix of Eq. 7 equal to 0 one obtains for the characteristic equation 


- ['^n ^22 ^22 1 - *^22 ^11 ^11 *^12 ^21 ^22 


( 8 ) 


Q 


1 - n 


“2 ^ 

K21 S ^2 <^11 r + T'TIJ ^22 '^ll ^^22 ^^11 ^22 ‘ ^12 ^ 21 ^ ^ ° 


since <22 " >^12 *^21 ~ ^ ' 

It is very tempting to drop cross-coupling terms like <12 S21 022* 

K21 S^2 8-|2 S21 in comparison with other terms, in which case 

Eq. 3 factors and the waves types are comoletely decoupled. However, 
not all of these terms are negligible in any given situation. It is still 
instructive to establish a few criteria that facilitate approximate solutions 

of Eq. 8 and identify the propagation constant for sound waves in the 

air r, and sound waves in the structure Tc when possible. For that purpose, 

one need only consider wave propagation in the positive x - direction. 

First, note that the inequality 




'^12 .^21 ^^22 1 - n '^21 ^12 ^11 ^ ^11 ^22 ^22 '^22 ^11 ^11 


( 9 ) 


is valid provided 


(l - Q 

I n 


<2 - * '> 


< < 
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where ^ei ^ Qa = , and Qgs 

convenient parameters for analysis. 

This inequality is valid for all highly porous materials regardless of 
the value of the elastic ratio and at all frequencies of interest 

in the audio range. Moreover, since for ks. < 4 and 

reasonable values for pg , frequency is no limitation. In the experiments 
to be described later the value of the left hand side of the inequality is 
about 0.007 for ^ ^ lowest value for was 2.5. 

Unfortunately, there is no precise information about fiiCs/Kg for the materials 
studied but the best judgement is that > 1 over most of the 

frequency range of interest. The conclusion is that all of the cross-coupling 
terms in the coefficient for in Eq. 8 can be neglected. 

The other requirement to be satisfied in order to eliminate cross-coupling 


entirely is that 


$12 $ 2 ^ < < $22 


( 11 ) 


in Eq. 8. 


This inequality reduces for 


')e1 ^ ’ 


to 


1 < < (Qef + 


( 12 ) 


which is more difficult to meet than Eq. 10 by orders of magnitude at 
low frequencies, say below 300 Hz. The requirement is met in the experimental 
situation at frequencies above say 1200 Hz. The conclusion is that the 
coupling terms 3^2 ^21 ^ safely .be neglected at all 

frequencies of interest. 



In order to obtain a solution to Eq. 8 valid over the entire range 


of frequencies of interest consider the requirement that 
4 ^22 ^11 <^22 ^^11 ^22 " ^12 ^21^ ^[*^11 ^22 ^22 ^22 ^ll^ll] 


Again, under the assumptions that < <1 and k-ii 1 this inequality 
is valid provided both 




(14b) 


which can be met under certain conditions. There are two extreme conditions 
in the relative values of the ratios '^es'^'^a satisfy 

both requirements. The first is the so-called "flexible blanket" frame 
condition where 


Qk. 


1 > > 2 Pg/Pgg 


(15) 


or Cg^/Cg^ > > 2 


where c 


2 _ 


(Pj,K,a) 


-1 


9 2 1 

and c (p K )" . Here c, is identified as the speed of sound in 
s '■^es s a 

the air of the pores and c^ is the speed of compressional elastic waves in 
the frame, Eq. 15 is easily satisfied for most highly porous flexible 
blanket materials over wide ranges of frequency. 
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The second extreme condition is the so-called "stiff frame" condition 

where 



> > Pes/2Pa 


( 16 ) 


or < < 2 which probably is a more difficult condition to 

satisfy over a wide range of frequencies. From equations 2 and 4, 

Pes ^ ^ (1 - n) ^ *^s ~ ^ , which for many porous materials including 

Pa Po 

flexible blankets is a number much greater than unity. For one type of Kevlar 29, 
the ratio p /p =40 based on a value of kg = 2.1 while for Scottfoam 

0S 3 

p Jp =13 based on a value of k = 2.1 and numerical values of 

Pg for each material (see Table I). It appears on the basis of measurements 
by others, notably Beranek^^^, that commercial materials which are acoustically 
stiff at low frequencies lose that quality at high frequencies, i.e. the 
ratio Qk /k is highly frequency dependent for some materials notably 

S d 

stiff acoustical tiles or boards. 

Solutions of Eq. 8 are quite different depending rather critically 
on choices for the stiffness ratio density ratio 

p /p . Again note that k-., contains ^^Kg/ic^ as a parameter 
but unless extremely large say greater than 50 and is 

adiabatic the value of i<ii will remain close to unity. If remains 

isothermal at all frequencies its value is unity. 

It is not worthwhile speculating further about the relative values of 
ic and K until more information about k for spun fiber and ooen-cell foam 

S3 ^ 

materials is available. It is doubtful that static values of will suffice and at 
best would provide a bound on * No concerted effort to find 

dynamic values of <g for highly porous, acoustical materials has been made 
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to-date. Also, such information is not provided by the manufacturers of 
acoustical materials on a routine basis. Instead, it will be assumed 
that the conditions of Eq. 13 obtain and then proceed to find 
approximate solutions for r . 

B. Propagation Parameters 

If the conditions of Eq. 13 are satisfied then the two solutions 
to Eq. 8 can be expressed as 


■^22 1 - R 


‘^11 ^22 


($11 $22 - Si2 $2l) 


Kii 022 ^22 




(17a) 


1 - Q 


^22 ^11 


*11 


and 


r 


2 

+ 


<11 a ^2 5^2 + 1^22 


1 - n ^11 



(17b) 


At this point one cannot identify whether r_ is associated primarily 
with wave motion in the elastic frame or sound waves in the air pores and 
likewise for r+ , However, if the conditions of either Eqs. 15 or 16 
are imposed then the identification is readily made by comparison with 
exact solutions for decoupled motion. 

It is convenient to rewrite Eq. 17a in the form 


0) 




c2 


(1 ^ r > " r > 

^es ^a 


^ei 

Qa^es 


(1 + 


Q • 

^ei 


(i " r ) 






( 18 ) 


Also, Eq. 17b becomes 


2 


r 


+ 
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2 2 

If conditions of Eq. 15 are satisfied then > > Cg and 



Note, in the completely decoupled case that the sound waves in the 
.pore space have a propagation constant 



and the elastic waves in the structure have 


( 20 ) 


(21) 


( 22 ) 



2 2 

Thus, for the flexible blanket where ^ ^ S * ^a 

identified with r_ and is identified with . 

The most important effects of wave coupling are identified by examining 

equations 20 and 21 and noting that the effective wave speeds and internal 

damping are modified by terms that involve the product QggQg* (See Appendix A 

2 2 

for a more precise analysis.) Since » Cg and K-j-j - 1 , the last two 
terms in Sq. 21 usually are dropped and the expression for Tg simpTifies 
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to the decoupled case. Also, Eq; 20, can be expressed in terms of Q and 


Pes/Pa 


as 



which is a convenient form for parametic analysis and numerical computations 
if coupled motion is important. 

All of the parameters in Eq. 24 are presumed known or can be otained 

from independent measurements. Values of $ for Kevlar 29 and Scottfoam 

are available^^’^’^^ . However, values of k are more difficult to establish 

s 

and Cg, Pgg/Pg and are dependent on that value. It is a very 

important parameter to establish in the theory of Zwikker and Kosten. No 
assumption about the value of has been made except that it has been 
treated as a real number. It is pure real at either very low or very high 
frequencies but at intermediate values it is a complex number. In the 
parametric analysis it will be presumed known at high frequencies, either 
adiabatic or isothermal. It is expected that if there is appreciable mechanical 
damping in the fibers of the frame that also will be a complex number. 
Fortunately, in the parameteric analysis a precise value did not have to be 
established. 

The "stiff frame" conditions of equation 16 is interesting and useful for 

better understanding of porous acoustical materials. If > > 2 p „/p,, 

as es a 

2 2 

say greater than 50, then c_ >•> c. and this theory predicts that 

b a 
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and 


c/ 


2 2 

r = r 

a 




(1 + 


) - ^ 

> 2 


(l-K,,) (l.i) 

where now 


'es 


^11 ^ ^ 


■'ll *>65 '’a ^ or > 

a 


(26) 


unless *^3 • Again, the express! 


on 


for r can be simplified by dropping the last two terms which is fairly 


good approximation if 


> > 2p^ /p^ 


JiKg ~"es' ''a 


Then T^ reduces to that 


for decoupled motion with slight modifications in that ‘^n ^ ^ 

Note, also that the roles of r_ and have interchanged from 
that of the "flexible" porous blanket. In effect, there is a switching 
regime due to internal wave coupling which is most likely a frequency dependent 
phenomenon. Just how and where this switching occurs depends on the detailed 
behavior of 

o a 

On the basis of experimental information now available it appears that 
<a 

the inequality ^ ^^es^^^a obtains for rigid board type materials at low 
s 2 


frequenci es i . e . Q, 


es 


< < 1 


At high frequencies it does not seem to matter 


whether the porous material is rigid or not and the waves decouple in any event. 
Predictions are then fairly readily made and there is good agreement among 
various workers. Delany and Bagley^^^^ have noted this and characterized a 
wide range of fibrous absorbent materials on the basis of a single parameter 
such as w/<J> . They used various grades of glass fibers and mineral wool 
materials and found an excellent empirical fit over a range of frequencies 
from 250 to 4000 Hz which includes the range used in the this study. However, 
at sufficiently low frequencies there is appreciable wave interaction and it 
is doubtful that such a simple characterization is possible. This, is undoubtedly 
the case where the structure factor kg is greater than 4. 
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( 2 ) 

Zwikker and Kosten' have discussed some of these matters in detail and 
conclude that the frequency above which coupling forces are likely to be in- 
significant is given by 



where is designated as a decoupling frequency. For Kevlar 29 using k^ == 3 

this frequency is about 130 Hz which is well below the lowest frequency used, 

namely 300 Hz. For Scottfoam it is about 85 Hz for k^ 2 

2 

It should be noted that for conditions where > > 1 , both the 
"stiff frame" and the "flexible blanket" approximations yield identical results, 
i.e., Eq. 26 becomes identical with Eq. 20. However, at lower frequencies, 

say Qgg < 3 there is considerable discrepancy, see Figures 1 and 2. 
Experimental facilities were not available to check the low range of 
frequencies (below 300 Hz) directly. for either Kevlar or Scottfoam. 

Shown in Figure 1 is a plot of the normalized propagation constant 
c /tor = B + i A as a function of Q based on Eq. 24 for 

p^^/Pa = 38 . Also shown for comparison purposes are results predicted 
byEq. 22 based on completely decoupled motion. Note,' the large 

discrepancies for the range < 4 especially in the attenuation component. 

On this basis any analysis based on decoupled wave motion will overpredict the 
attenuation and underpredict the phase shift primarily at low frequencies. It 
also should be noted that the greater the ratio Pgs/Pg the more sensitive the 
motion is to internal wave coupling. Further, it should be noted that the real 
and imaginary part of at a given frequency are not in general equal unless 
the waves are decoupled and is low, say less than 4. In particular, 

A and B are not equal for large Q . This judgement is based on the 
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assumption that o-./p.. > 10 . 

cb u 

In closing this discussion it is interesting and useful to point 

( 4 ) 

out that the general results are compatible with earlier work by Beranek 
except that he employed some different definitions. He studied the 
properties of both "flexible blanket" and "rigid tile" and expressed his 
results in terms of an "effective density" which is different from 
p and a dynamic resistance here called . This matter is discussed in 

the appendix. The two analyses are brought into better Quantitative 
agreement if p^Jl becomes ^p-|^ in Beranek 's work. The qualitative results 
gy >0 for all intents and purposes the same. 

In the following section the acoustical impedance and sound absorption 
coefficient relationships for bulk materials will be discussed briefly. 
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III. Impedance and Absorption Relationships 

The utilization of the coupled equations of motion to predict acoustical 

( 2 ) 

impedance of absorptive materials has been discussed by Zwikker and Kosten ■ , 
Kosten and Janssen^^^^ and Beranek^^^ in great detail. Generally speaking 
the analyses are very complex and difficult to follow in detail. Here, an 
approximate method is employed which yields good results in the "weak 
coupling" cases. In such cases the coupling terms in the dynamical equations 
of motion are dropped completely and the impedance for the sound waves in 
the air of the pores becomes 


Z 


a 


wp. + i$ 

u 


(28) 


on the basis of Eq. 6a and for waves in the elastic frame 


Q 

-s "" (1 - flT 


toPes + 


(29) 


on the basis of Eq. 5a. 

Any effects of wave coupling then enter only through the modified values 
of r and r 

a S 

According to Kosten and Janssen^ this can be done quite rigorously 
when the frame is either very flexible or very stiff by very careful 
approximations. But, apparently this procedure cannot be followed in 
intermediate situations where c, = c^ . Thus, if the condition imposed 

a S 

9 ? 9 9 

is either c >>c ore <<c_ then Eqs. 28 and 29 can be 

3 S 3 S 

employed to determine the impedances. 

Here it will be assumed that the "flexible blanket" condition obtains, 

Eq. 15, and one predicts that 
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The form of Eq. 31 is not particularly useful since we do not 

2 2 2 , 2 

have information about the ratio /c^ . So usually /c^ is 

taken small compared with unity for flexible blankets and then for reasonable 
values of Q and Q the last two terms in the denominator of Eq. 31 are 

G$ 3 

dropped. 

The total surface admittance of a sample of bulk porous material is 
Y = Y + Y, where the area weighting has been taken into account. 

Since o^c > > o^c one can establish that |Y 1 < < [Y | and 

the admittance of the sample is that due to the air motion alone. 

Obviously there are restrictions on both "S^^a 

P Jp for this to be the case. The approximation presents a great 

GS 8 

simplification and it usually is made in most applications of porous 
materials and in other applications like predicting the influence of 
ground on the absorption and reflection of sound. 

Then the normalized impedance of the surface of a bulk material becomes 
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\ X ~^) (32) 

a [l +/te\ — 1_ + i/P^\ ^es 1^/2 

I \ Pa ' Qes + 1 \ ^a/ Qes^ + 1 J 

where is the adiabatic speed of sound and the parameter k = k^/Kg^QY • 

First, one notes that in general quite a large number of parameters are 


required to characterize the impedance properties of the porous material. 
At least four standout explicitly, namely k, SI, Pgg/Pg’ ^es ' 

Recall that p^3/p^ = "^p^ \ f - ^ • Also = «> <>es/^ 

requiring knowledge of the specific flow resistance $ . This is a larger 


number of parameters than usually employed. 

The structural factor kg enters into all of the parameters namely 
k, Pes/ and However, if wave coupling is neglected then Eq. 32 
simplifies and can be expressed as 



?a_.,Xa_, VL(l.iJ_)’^^ 

Pc^c.^ P.^c,^ Qa 

0 0 0 0 “ 


(33) 


where = aipQ(ks/Sl<J>). 

Now kg enters into only two parameters, k and Q^. The high frequency 
asymptote in both Eqs. 32 and 33 is Vk/fi • I"*" ^a known 

at high frequencies then k^ can be established uniquely, in principle. 

The normal incidence sound absorption coefficient can be expressed 


in terms of the surface impedance and becomes for highly porous bulk samples 
namely 

A (Rg/p„c„) 

“a= — . (34) 

(1 *Wo)^ ^ (VPnCj^ 


a' '^0 0 ' 


Generally speaking 0(3 requires as many parameters as in order to 

" a 


describe it behavior. 
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If normal impedance measurements are made on samples of porous material 
of finite depth backed by a plate that presents an acoustic impedance 
large compared to P^Cg then the impedance of the surface is modified 

to read 


iuL 

f’o‘^0 


bo^o 


' TO _ _a_ 
'Vo ' 


P C 
^0 0 


coth^ir il) 

a 


(35) 


where fa ' ^ “a e°"'PeWt1onal purposes Eq. 35 can be 

expressed as 


Rin (f^a^‘^o‘^0^ coshi'a^J,) + (Xg/OgC^) sin(6g£) cos(6g^-) 


p c 
*^0 0 


cosh^ - cos^ 


(36a) 


and 


Xin (^a'^^o^o^ ^^a^^ cos(3gJi) 


P c 
^0 0 


cosh^ (otai.) - cos^ (^a^^ 


(36b) 


If, at a given frequency *^a'^ taken so large that tanh (cx^^-) ^ 0.998, 

which obtains for a^5- i 3.6 , then for all practical purposes 

R = R. and X = X. .This condition was met or exceeded for the 
3 in 3 in 

range of frequencies used in this study. 

A representative plot of f^a^^oS Kevlar 29 

is shown in Figures 2a and 2b, respectively. Again note, the range of 
for which the decoupled wave approximation is valid. 


In the following section a method for establishing the value for the 
structure constant k^ is presented. Also, a brief summary of some previous 
work endeavoring to relate k^ and will be given. 



IV. Parametric Studies 


While there are a number of ways of describing the impedance and 
absorptive properties of bulk porous material one would prefer to employ 
as few parameters as possible. In particular it is advantageous to 
characterize both the real and imaginary part of the impedance so as to 
be able to predict just what physical properties essentially control their 
behavior. One then could select optimum values for all parameters in 
order to achieve a desired impedance or absorptive characteristic, if possibl 
Also there are hidden parameters such as temperature to which the impedance 
and absorption could be sensitive. 

Some attempts have been made to characterize the impedance properties 
of fibrous absorbent materials by a single parametric ratio such as w/$ . 
This may be possible if the structural characteristics of the candidate 
materials all have about the same values. 

The impedance of a bulk sample of a porous flexible material is given 

in Eqs. 32 and 33 where the high frequency resistance is 



In the theory of Zwikker and Kosten Pa^^o ~ "^a is assumed 

known at high frequencies. 

Some workers prefer to let Pa'^Po " ^ where m also is interpreted 


as a structure factor. 
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where k_ = m if k, = — — and k - “ if S ” 

^ ^ POY a Y a *^0 

characterizations m is interpreted as a function of the porosity alone. 

Note also that the speed of sound in the air of the pores 

c 2 = _J — = c ^/k 
*^a K- p.n 0 'S 

a a 

2 _ 1 

where is the speed of sound in the air of the pores if the 

air motion is unrestricted. However, there appears to be no immediate 

advantage in characterizing behavior of the material by simply redefining 

o /o so the definition in Eq. 2 is retained. 

It is surmised that the parameters $ , ^ and k^ are interrelated 

but the question as to how probably cannot be settled as a purely empirical 

basis. In the theory of Zwikker and Kosten, the specific flow resistance is 

independent of frequency in the audio range. A precise definition of low 

and high frequency behavior is not available for most materials but it 

appears for highly porous, fibrous materials including open cell foams that 

the acoustic Reynolds number wd /v must be less than unity in order for ^ 

to be independent of frequency. Here d is a mean diameter of the fibrous 

filament and ,v is the coefficient of kinematic viscosity for air surrounding the 

fibers. The average spacing between fibers must be sufficiently great so 

that there is no interaction between fibers. If the open-cell foam consists 

of solid filaments of mean diameter d with a mean separation D»d, then 

M3l 

It is likely to have a comparable index' . 

On the basis of experimental data available it appears that 0 for Kevlar 29 

and Scottfoam is constant in the frequency range below 3000 Hz. However, at 

2 

high frequencies, where tod /v > 1, it appears on a theoretical basis employing 

simple acoustic boundary layer theory that <1> should vary as v/ to . 

(5) 

R. W. Morse' ' has predicted for acoustic wave propagation in porous media 
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modelled as ensembles of small tubes that 


IT- “ ^Po 


/d^, lod^/v < < 1 


and 


(40a) 


l/Zpp u)^ 


(jjd /v > > 1 


(40b) 


where d is the diameter of the tubes and ^ ” Pq ^ 

The constant of proportionality in Eq. 40 is important and likely to 
vary from one material to another. More recently Hersh and Walker^^^ have 
shown in the basis of a one-dimensional layered fiber model for Kevlar 29 


that 


I = (vpp/d^) 

S 



(41) 


where K = .059 is an empirical constant. In their work it does not appear 
possible to isolate "J* , , and k^ . This representation then agrees fairly 

well with their data over the frequency range 100 to 2000 Hz. It is interesting 
to note that the theory of Zwikker and Kosten can be brought into substantial 
agreement with the work of others by such an identification but generally 


there is an empirical constant to fit. In the parametric work of this 
paper values for $ and n were considered known and a single parameter Icg 

selected to provide a fit of the data over the entire frequency range. 
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In order to establish a value for kg parametric plots of Eq. 33 based 
on decoupled motion were made employing as a variable and y/k/Q as 
a parameter. Initially was selected to be adiabatic and 
determined from the high frequency asymptote. That provided an initial 
estimate of k . Then Q, values were calculated and the data plotted 
against the parametric plots. A new value for y/k/^ was then established 

on the basis of comparisons of the ^a'^^oS “a 

with the predictions. This provided a correction to the k^ value. It was 
not possible to choose a single value for k^ that fit both the high frequency 
and low frequency data for all types of materials. One reason may be that Xg 
varies from its isothermal to adiabatic value as frequency increases. The 
most satisfactory fit to the low Kevlar and Scottfoam data, especially 
the reactive component of impedance, was obtained by selecting xg to be 
isothermal at all frequencies and then slightly adjusting k^ to fit the low 
frequency data. This procedure also is the most satisfactory for fitting 
sound absorption data. Fig. 12. The fit to the resistive component at all 
frequencies is then not so good as desired. Another reason may be wave coupling 
at the lower frequencies that invalidates the parametric procedure. There 
are some noticeable differences between Kevlar and Scottfoam in this respect. 

No attempt was made to alter the values of ^ and in the procedure, only k^ . 

The parametric- curves based on Eq. 33 are shown in Figs. 3, 4, and 5. The 
range of parameters 0.9 <Vk/il £ 2 was preselected because most highly 
porous fibrous materials would fall in this range. Some fibrous materials 
having large values of k^ and low values of n may fall outside but are 
not likely candidates for most noise control applications. It may be possible 
to develop materials for which Vk/R < 1 which is highly desirable at 
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low frequencies, but they would require low values of and high values of i. . 
It is doubtful, in any event, that VlT/il would ever be less than i/y^ = *845 
for air. Also if one were interested in predicting the absorptive 
properties of porous ground cover perhaps a different set of Vk/n values 
would have to be considered. 

Finally, parametric curves based on coupled wave motion could be 
plotted against Q _ in other applications especially these involving extremely 
low frequencies, say Q . ^2 or 3. Such parametric curves might be 
more useful for engineering design where specific values for either 
R /d c or X /p c (or both) are desired from a bulk absorbent 
material. No attempt has been made to fit the measured data to a parametric 
set of curves based on Eq. (32) as of this data. Parametric curves should 
be developed for coupled, Eq. (32), as well as decoupled, Eq. (33), wave motion 
Results based on the parametric studies are presented in Table I. 

In the following sections the experimental results will be presented 
together with a discussion based on comparisons with decoupled wave motion. 
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V. Experimental Results and Discussions 

An experimental program was undertaken in order to determine the 
impedance and absorptive properties of two types of Kevlar 29 that were 
available. One type is a light weight blanket that appears to have a 
layered structure that can be easily separated. The other type is a more 
dense blanket that appears to have a woven structure and cannot be separated. 

The fibers appear to be the same in both types. The acoustical properties 
of the two types are quite different, the light weight blanket has a 
comparatively low value of flow resistance while the dense blanket has 
a higher value of see Table I. Values of $ for Kevlar had been 
established earlier by Smith and Parrott^®^ and Hersh and Walker^^^. 

Also available were samples of an open-cell foam material called 

Scottfoam. Its propagation and impedance properties had been studied and 

( 8 ) 

were available for comparison purposes' . No detailed comparisons and 
discussion of propagation contains for either Kevlar or Scottfoam will be 
presented in this report. It did turn out that the theory of Zwikker 
and Kosten overpredicts the attenuation and underpredicts the phase constant 
slightly in the frequency range 300 to 3000 Hz for high $ Kevlar 29 based on the 
parameter values presented in Table I. Experimental values of attenuation 
were used to determine the required depth of sample based on Eq. 36. The 
required depth of sample for Kevlar 29 was 12 cm while that used for Scottfoam 
was 25 cm. 

No detailed description of the experimental apparatus and techniques 
employed will be presented here. All impedance and absorption measurements 
were made in a conventional impedance tube having an inside diameter of 5.72 cm. 
The samples were held securely in a special tube andfaced with a highly 
porous screen at one end and a rigid impervious plug at the other. The 
technique employed plane waves to measure the standing wave ratio and positions 
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of the first minimum in the acoustic pressure standing wave patterns. The 
sound pressure level at the face of the sample was detected by a reference 
microphone and held constant. The standing wave detector was controlled 
by a dedicated digital computer which greatly facilitated the gathering of 
data and subsequent data reduction. 

The measured values of normalized acoustical impedance and normal 
incidence absorption coefficient are shown in Figs. 6-12. The trends and 
levels are clearly established for the three specimens tested. The data 
were highly repeatable and followed typical patterns for porous, flexible 
blanket type materials. Some difficulties in holding the dense Kevlar discs 
were encountered but were overcome by compressing the samples slightly. 

This did have the effect of raising the high frequency asymptote '/k/n from 
about 1.5 to 1.9 and corresponding reducing the absorption coefficients 
at all frequencies. 

The parametric study established the values for the structure factor 
kg shown in Table I. Other physical parameters used in the computations 
are included in the table for comparison purposes. The value of ^ for the 
light weight blanket Kevlar was computed from the data of Hersh and Walker, 
namely 46.4 cgs rayls. Their force equation neglects inertial wave coupling 
entirely by setting k^=l. They then proceed to develop Eq. 41 where their 
structural parameter K = 0.059 . If one then identifies $J^/ks with their 
parameter a the computed value is obtained. It is believed that the two 
formulations are self-cosistent if this identification is made. When preparing 
samples of the low $ Kevlar the porosity was preselected as 0.94 to match 
the high $ Kevlar and then the required amount of material was packed into the 
.specimen holder' to a dtepth of 12 cm. 
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The Scottfoam sample is an open-cell foam that under comparative 
examination with Kevlar has a fiber diameter about 1.667 times as great. 

If the published value of d = 1.3 x 10"^ cm for Kevlar 29 is correct then 
d for the Scottfoam is about 2.1 x 10"^ cm. This value for d then 
correlates very well with a measured value of $ = 18 cgs rayls. 

The characteristics of Rg/Po^o’ °^a computed 

using Eg. 33 and the values of the parameters in Table I. No further 
adjustments in any of the parameters were made and the results of these 
computations are plotted together with experimental data in Figs. 6-12. The 
predictions of X /p c agree very well with measured values for all 
samples. There are, however, some discrepancies in the predicted values of 
Rq/PqC^ especially at low frequences. The exact reason for this 
discrepancy which seems to be outside the range of experimental error is 
unknown as of this date. However, it was decided not to adjust k^ to provide 
a better fit to R,/p„c^ until a parametric study based on coupled wave 
motion is undertaken and completed. 

Predicted values of absorption coefficient ot^ are sensitive to small 
changes in X^/p c„ and provide another index for comparison purposes. The 
agreement between theoretical values and experiment is excellent for both 
the low $ Kevlar and Scottfoam over a wide range of frequencies. The low 
$ Kevlar characteristics do depart in an intermediate frequency range providing 
some evidence that the samples were vibrating and absorbing energy. 

Of special interest are the impedance and absorption characteristics of 
the high $ Kevlar, especially the larger value of k^ required to 
provide satisfactory agreement. This result was not entirely unexpected 
since the value of $ is higher also. Note, that <I>.and k^ are directly 
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f 5^ 

proportional on the basis of theory of Morse, Eq. 40. Moreover, Morse' ’ 
predicted that the structure factor kg will have a value of about 3 if the 
fibers are randomly oriented. He also reported values for kg between 2 and 3.4 
on the basis of measurements in granular porous media. 

The other conclusion based on the parametric studies was that the thermo- 
dynamic process is essentially adiabatic at high frequencies for the high 
Kevlar. Selecting an isothermal value for , while providing a better 
fit to the R,/p c data, did not agree well with the X,/p c^ data 
being much too high, and correspondingly produced values of that 
were much too low especially at high frequencies. It is concluded that the 
thermal process changed with frequency and a frequency dependent value for 
Kg is required to provide a better fit of the data. It can be remarked 
with confidence that the two types of Kevlar behave quite differently both 
with regard to impedance and absorption. This matter would require more 
theoretical study based on very careful thermodynamic and structural 
considerations to provide a better rationale for the differences noted. Note, 
that both types had the same porosity. 

On the other hand, the low Kevlar and Scottfoam samples do behave 
surprising alike in spite of quite different mechanical structure. The 
Scottfoam has a much lower $, a lower material density , and a lower 
structural to air density ratio Pes^^a* However, the Qg of Scottfoam is 
roughly twice that of Kevlar for the same value of k^ . This fact makes it 
a more efficient absorber at low frequencies even though is the 

same. In Fig. 12 the frequency range for Scottfoam is 500 to 3500 Hz while 
that for Kevlar is 300 to 3000 Hz. 

When plotted versus the data of low $ Kevlar and Scottfoam 
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collapse to the same parametric curve. This is shown quite convincingly 
in the absorption characteristics. Fig. 12. The high Kevlar data, on 
the other hand, follow a higher parameteric curve and display a 

correspondingly lower a characteristic. The ratio is justl.9 to 1.3. All 

d 

of the absorption data could be collapsed when normalized with respect 
to this ratio. 

It is becoming common practice to plot absorption data versus a 
parameter such as The values of are then varied by 

changing frequency since ^ regarded as constants. The 

parameter can be interpreted as an energy ratio, namely of 2ir times the 
maximum stored acoustic energy per unit volume of material to the energy 
dissipated per unit of time per unit of volume. Thus, it is very convenient 
parameter from an energy viewpoint. It also can be expressed as 



where p c is the acoustic impedance of the sample neglecting dissipation 
and Aa is the acoustic wavelength within the material. Other interpretations 
may be possible. Therefore, it was decided to plot all data and make 
comparisons employing as a frequency variable. 
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VI, Suninary and Recorronendations 

A, Results 

(1) A theoretical formulation of acoustic wave propagation in porous 
elastic structures due to Zwikker and Kosten has been recast into a form 
that makes identification of various coupling parameters notably inertial, 
frictional, elastic, and dilatational more evident, see Eq. 7. 

(2) Conditions under which the compressional elastic wave in the frame 

decouples from the acoustic wave in the air of the fibrous pores have been 
carefully examined. Weak coupling solutions have been formulated and 
expressed in terms of a frame structure quality factor which provides 

a convenient index, see Figs. 1 and 2. 

(3) A density ratio p /p has been defined which plays an important 

GS di 

role in the theory of wave coupling. The lower the values of this ratio the 
less susceptable the air acoustic wave is to the influence of energy storage 
in the elastic wave and vice versa. This ratio in conjunction with the stiffness 
ratio fiK /<, controls the nature of the wave coupling, see Eqs. 15 and 16. 

(4) A new derivation of the expressions for equivalent density Pg and 
dynamic resistance for porous materials has been developed. An equivalent 
quality factor also has been defined. This formulation could be used for a 
parametric analysis of materials where wave coupling is involved (see Appendix A.) 

{ 5 ) Parametric curves for impedance and absorption coefficient based on 
decoupled wave motion have been prepared which facilitate identification of 
the structure factor k^ for flexible bulk porous materials, see Figs. 3, 4, and 5 
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(6) Experimental characteristics of normal impedance and absorption 
coefficient in the frequency range 300 to 3000 Hz have been established 
for two types of Kevlar 29 porous blanket, see Table I. 

(7) Acoustic parameters based on the theory of Zwikker and Kosten for 
Kevlar and Scottfoam have been established, see Table I. 

I 

s 

(8) Agreement between measured and predicted values of impedance arid 
absorption ceofficient based on bulk properties is satisfactory. The wave 
motion was assumed to be decoupled but it appears that some coupling 
effects are in evidence, see Figs. 6-12. 

B. Recommendations 

(1) No concerted effort to refine the characterization of porous, flexible 
blankets for critical applications has been made. That might involve preparing 
a set of parametric curves and procedures based on coupled wave motion. If 
Kevlar or Scottfoam or any other porous, flexible blanket type material becomes 
a candidate material for duct liners then such an effort may be warranted. This 
also may be an essential procedure if effects of extended wave reaction on 
impedance of porous, bulk duct liners is tobe understood and predicted. 

(2) If angle of incidence effects on impedance and absorption of duct 
liners or other porous materials are to be fully predicted then a comprehensive 
theory based on a three-dimensional wave model is required. Once the essential 
parameters have been established then wave coupling can be introduced and its 
effect on the surface characteristics predicted probably via numerical analysis. 
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(3) Effects of finite sample size are not included in the theory of 

Zwikker and Kosten as applied to the testing of bulk materials in an impedance 

tube. Therefore, measured values of impedance and absorption can differ in 

detail from predicted values even if the acoustic parameters are essentially 

(15) 

correct. Such discrepancies have been noted in work by others ’ and 
similar discrepancies also were found here for some samples. For example, 
if the samples of high $ Kevlar are not mounted firmly in the holder and 
held by compressional forces they will vibrate and resonate at particular 
frequencies by virtue of their finite diameter. This is evidenced by extreme 
changes in impedance over fairly narrow frequency ranges which are not 
predicted by theory in detail. Even with precautions, samples of Kevlar 
exhibited some irregular behavior as noted in experimental values, see Figs. 8, 

9, and 12. 

(4) The elastic bulk modulus must be carefully defined and measureable 
if accurate information about wave speeds, wave coupling action, and extended 
reaction effects in porous materials is required. Good definitions, methods of 
prediction and measurements of all parameters is a goal worthy of further 
study. 

(5) The theory of Zwikker and Kosten can be used to characterize the 
acoustical impedance of porous ground. If characterized properly such 
information could be used to predict "excess" attenuation of sound attributed 
to absorption and reflection in applications such as aircraft fly-over noise, 
freeway traffic noise, the certification of noise sources, industrial plant 
site noise, and a host of other community noise prediction methods. This 
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effort might require both an extensive experimental as well as theoretical 
study to establish the appropriate characterization. 
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APPENDIX A 

CONCEPTS OF EQUIVALENT DENSITY AND DYNAMIC RESISTANCE FOR POROUS MATERIALS 

The purpose of this work is to derive expressions for the equivalent 
density Pp and the dynamic resistance $ applicable to flexible, highly porous 
blankets and rigid, highly porous tile. These are extreme cases that can be 
handled on the basis of the analysis of Zwikker and Kosten. 

The proper procedure is to examine expressions for the wave propagation 
constant r in cases v/here wave coupling is significant and then endeavor to 
express it in the form 


(l + i ^ ) 
\ wpp y 


A-1 




which is valid for decoupled motion where 

In order to do this in a systemmatic way one invokes the inequality 


'12 *"21 22 1 - 21 12 11 


" ^ 11^^22 
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A-2 


Which easily holds for highly porous materials, i.e. n 0.9, having a structure 
constant k^ _< 4 for all frequencies of interest, see Eq. 13. Also, 
the value of the stiffness ratio is not critical here so that A-2 

is valid for both extremes namely /k, » 1 and fiK./K:,.i.<< 1 . 

S a S a 

By writing the expression for r , Eq. 11, in the form 


r ^ - b + c = 0 
e e 

one finds after invoking A-2 that 


b 



+ 


Ksa 


sa 



^11 + Ksa 


sa 



A-3 


A-4 


where 


and 


Also one finds that 
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where the Q's have been defined in the text. 

2 

On the assumption that Qg^ < <1 at all frequencies of interest the 
ratio c/b becomes 
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where 


^sa " '^sa^/Qes 


and 


a - f ‘^SQ- VTr^ 

< -"sa ?%s 

Note, if kg = 1 > then = 0 and there is no inertial coupling 
For flexible, porous blankets, t.. << 1 and 


sa 

^sa ^ ^11 ^ ^ \ Q ~ 
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which reduces to 


(^) <’ -sV ' 


an inequality that is almost identical with the inequality in Eq. 13 and 
hence easily satisfied. Thus, from A-1 
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where 


and 
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(1 + i 


e' ^e 
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since Pa^^es ^ ^ ’ 

Pa 2 

At high frequencies, Q . » 1 and p = p while at low 

2 ^es ^ ■ 

frequencies, i.e < <1> the effective density of air in the pores 

of the material is p^ + p^^ indicating that the motion of the air 
particles and the fibers are tightly coupled together and move as one 
entity. 

It also is interesting to note that an equivalent quality factor 
can be defined as 


Q = = Q (1 + _es 

^ P 

e 


1 


) 


A-9 


^es 


which can be approximated by only on the added assumption that 


"es 


> > p /p, 

^es ^a 


This implies that any empirical analysis based upon a 


single descriptor like Qg must be justified carefully and likely to have 
limited applicability especially with regard to frequency range, most likely 
high frequencies. 
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2 

At low frequencies <J>g = <J> for flexible, porous blankets . 

indicating that the dynamic resistance is very low and varies with the 

2 

square of the frequency while at high frequencies, where Qgg > >1, 

becomes identical with the specific air flow resistance $ since 

no energy dissipation in the fibers has been taken into account. This result 

does not imply that ^ is always independent of frequency, however. These 

( 4 ) 

formulas agree fairly well with Beranek's expressions' if n replaces 
the equivalent density in his work. 

The other extreme where obtains for rigid porous tile or rigid 

fibrous boards. In that case 
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which is of the desired form. One then readily identifies 



(^11 * "^sa] 


or 



A-11 


which is a ratio that is sometimes called a structure factor and probably 
fairly independent of frequency. This result was verified by Beranek for 
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wire screen mesh but not for commercial tiles due to their inherent 
inhomogeneity. However, the behavior of with frequency for rigid tile 
materials and rigid fibrous board is not known. If the air motion is 
isothermal (at very low frequencies) 


Pg = Pq 

for Q ~ 1 and kg. <_ 4 independent of the value of . Thus at low 
frequencies - p^ the apparent density of the air of the pores. Hence, the 
frequency behavior of Pg for flexible blankets and rigid board is quite 
different. 

The equivalent dynamic resistance for rigid board is 





) 


A-12 


There is no apparent frequency variation in for rigid board either. 

The second term nKg^'^a • ‘^'^11 at all frequencies. However, Beranek 

found that 4>g decreased with frequency at sufficiently high frequencies for 
reasons unknown. At low frequencies he showed experimentally that $ for 
fine wire mesh screen. This suggests that Eq. A-12 is a low frequency formula 
and hence not valid over a wide range of frequencies. It might be concluded that 
rigid screens can become acoustically flexible at high frequencies and are stiff 
only at low frequencies. This matter bears more careful study and perhaps a 
reformulation of the dynamical equations of motion is in order to account for 
material stiffness. One possibility based on a casual observation of cross 
coupling terms in Eq. 10 is that no dynamic stiffness cross coupling terms 
appear, i.e. o ■^2 = Ogi =0 in the coupling matrix. Perhaps such coupling 
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terms should be included in order to develop a more complete theory applicable 
to a wide variety of porous materials and valid over wide ranges of 
frequency. The theory as developed by Zwikker and Koster does appear, therefore, 
to be applicable to highly porous, flexible materials and rigid porous 
materials over a limited range of frequencies. Hence, it has limited 
applicability in so far as engineering design of materials is concerned. One 
would desire to be able to predict both the real and imaginary parts of 
and specify them in terms of material parameters at all frequencies of 
concern. 

In the theory of Zwikker and Kosten knowledge of the frequency behavior 
of Qk_/k^ appears essential and should not be limited to predictions of unper 
and lower bounds. This is an area that deserves more study. 

Finally, Eq. A-1 as here formulated could provide the basis for a more 
systemmatic parametric analysis of a wide range of porous materials including 
flexible blankets, ground cover, and stiff tile. This also is an area that 
deserves more study. For example, it is not known whether the inequality 
of A- 2 is valid for various types of porous ground but it appears to be easily 
satisfied for a wide variety of commercial materials designed as good absorbers. 
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APPENDIX B 


FUNDAMENTAL EQUATIONS OF MOTION 

( 2 ) 

In this appendix the fundamental equations of Zwikker and Kosten' ' 


will be recast into the form employed in the text of their report. In an 
elastic porous solid containing air the equations of motion in the elastic 
frame are given as 


and 


3 Pi 
9X 


= P, “—L + s 


!fi = K _ 

3 1 ^1 3 X 


3'>, 

1 -^-1- 2 

(1 -h ) 
h 9t 




and in the air of the pore space as 


^2 _ 


9 v*" 


and 


9x 

9 P 2 
9t 


“ "2 3“r'" = *''2 ■ '’P 

9 

h K? (^2 

at 


9 Vi 

P ) — - 
0 ^ 3x 


a) 

b) 

B-1 

c) 

d) 


where the parameters p-| i P 2 » i<i . '^ 2 ' ^ defined on page 52 

in Zwikker and Kosten i In their work p^ is the force per unit area of cross- 
section of the sample acting as the frame, Vj is the mean velocity of 
the frame, P 2 is the excess force per unit area of cross-section of 
the sample acting on the air, and V 2 is the mean velocity of the 
air in the frame. 

For convenience the Eqs. of B-1 are recast in operational notation as 


9Pi 


|[p^ + P2(l<s - 1 )]^+h2a}v^ - [p2(k^ - l)3-T+h^^]v2 


3 .u2 


a) 


B-2 


9P 
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Pi o 9 Vi 3 V, 

^ = [K^ + (1 - h)2 (,C2 - P,)] ^+ (1 - h) K2— 


b) 
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'2 ""s 9 t 
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c) 
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9 Vr 


9 Vi 


-Tt = ''*'2 IT <^2-V-TT 


d) 


after substituting for the coupling coefficient s. 

In order to write the parameters and variables in a more common 
notation define: 


h = Q 

p, = (1 - n) Pj 

P 2 = !2 Pj, 

^2 *" 1/^ 

Vi = 


a = 0 

P1 = (1 - n) Pg 

P2 = ^pQ 
*<2 = 1/^a 
K, = l/<3 


B-3 


Here p is interpreted as the dynamic pressure of the solid frame, p 
is the acoustic excess pressure in the air pores, is the mean velocity 
of the frame, and v is the particle velocity of the sound waves in the air 

d 

pores. 

By assuming harmonic excitatti on at an angular frequency to and propagation 
of elastic and sound waves in the x-di recti on according to the convention 
exp [-i(r X - wt) ] Eq. B-2 can be expressed as 


r Pg = n (o3 Pg 5 + i $) Vg - ( to pg^. + i 0) 

(1 - n) P3 = r V3 f .,2 r 

I Pg = (topg + i $) - Q (lopg. + i$ ) Vg 

“^a ^ ^a~ '^22 ^ ^ *^21 ^ '^s 


a) 


b) 


B-4 


c) 


d) 


and 
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where Pg» Pg> and are complex amplitudes and the parameters 

n n n IT e- v- and 9*^6 defined in the text. 

^ Pes’ Pei* Pa’ IV 12’ ^22’ 21 

It is worthwhile noting that Eqs. B-1 (b) and (d), the equations of 
continuity, have been criticized and reformulated on various grounds^^^^. 

Also, the definitions of p-j and P 2 render these equations invalid^^®^. 

The upshot of this matter is to simply remove the term P k from 

0 cl 

the dilatational coefficients and < 2 i in Eqs. 5b and 6b of the text. 

However, for highly porous frames considered in this report this procedure does not 
appear to be a significant change and it was decided to leave the definitions 
of Zwikker and Kosten intact. If the basic equations are reformulated for 


precision design purposes then such refinements may be warranted. 
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^\^terial 

Kevlar 29 

Kevlar 29 

Scottfoam 


low 

high ‘I’ 

open-cell 

Parameter 

Type 

Type 

Type 

<I> (cgs rayls) 

46.4 

79.4 

18 

*^s 

2.1 

3.2 

2.1 

n 

.940 

.938 

.950 


1 

1 


high frequency 


value 

’’o 

O 

□l- 

P 

0 

Ps (g/cm^) 

1.44 

1.44 

0.60 

^es'^^a 

38 

25 

13 

d(cm.} 

1.3 X 10"^ 

1.3 X lO'"'^ 

2.1 X 10“^ 


Table I. Mechanical and Acoustic Parameters of Test Samples. 
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Figure 1. 

Figure 2a. 

Figure 2b. 

Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 

Figure 7. 

Figure 8. 
Figure 9. 
Figure 10. 
Fi gure 1 1 . 
Figure 12. 


FIGURE LEGEND 

Plots of the normalized propagation constant-^ Tg = B +i A versus 

Qes "For Kevlar 29 based on Pps/Pa ^8. The subscript 1 refers 
to coupled wave motion and the subscript 2 refers to decoupled 
wave motion. 

Plot of the normalized resistance Rg/Po^ofor low Kevlar 29 
based on 38 and /k/JJ = 1.3. The subscript 1 refers 

to coupled wave® motion and the subscript 2 refers to decoupled 
moti on 

Plot of the normalized reac^nce -Xa/^CQ for low $ Kevlar 29 
based on Pes'^Pa " = l.T. The subscript 1 refers 

to couple wave motion and the subscript 2 referes to decoupled 
moti on . 

Plot of Rg/OgCo versus Qg employing Vk/n as a parameter for 
decoupled motion. 

Plot of "^a'^PoS versus Qg employing Yk/n as a parameter 
for decoupled motion. 

Plot of normal incidence sound absorption coefficient versus 

Qa employing as a parameter for decoupled motion. 

Plots of Ra/*^o^o versus Qg for low Kevlar based on parameters 
in Table I. Adiabatic bulk modulus is shown for high frequencies. 
Experimental data shown for frequency range from 300 to 3000 Hz. 

Plot of -Xa/PjCo versus Qa for low ^ Kevlar based on parameters 
in Table I. Experimental data shown for frequency range from 
300 to 3000 Hz. 

Plot of R,/p c versus Qa for high Kevlar based on data 
in Table I. ^ ° ° 

Plot of ~^JPo^n versus Qg for high $ Kevlar based on parameters 
in Table I. 

Plot of ^a^Po^o versus Qa for Scottfoam based on parameters 
in Table I. Experimental data for frequency range from 500 to 3500 Hz. 

Plot of "^a^Po^o versus Qa 'Fo'r' Scottfoam based on parameters in 
Table I. 

Normal incidence sound absorption characteristics for Scottfoam, low 
$ Kevlar 29, and high Kevlar 29. 
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Figure 2a.- Normalized acoustic resistance for Kevlar 29, r-i - coupled motion, r« - decoupled motion. 
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Figure 7.- Comparisons of the nomalized acoustic reactance for low Kevlar sample. 
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Figure 8.- Comparisons of the normalized acoustic resistance for high $ Kevlar sample 
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Figure 9.- Comparisons of the normalized acoustic reactance for high $ Kevlar sample. 
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Figure 10.- Comparisons of the normalized acoustic resistance for Scottfoam sample. 
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Figure ll.- Comparisons of the normali 
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Figure 12.- Comparison of the normal incidence sound absorptive coefficient of Scottfoam, low $ Kevlar 
and high $ Kevlar samples. 
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